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We calculate the quasiparticle scattering rate in the superconducting state of the overdoped 
cuprates, in the context of the Eliashberg formalism for a Fermi liquid with strong van Hove singu- 
larities close to the chemical potential. For a d^2_y2 superconducting gap, we demonstrate analyti- 
cally that the scattering rate is linear in the maximum of temperature or energy, but with different 

00 I intercepts and momentum dependence, thus extending our earlier results on the normal state. We 
discuss our results in view of angle-resolved photoemission experiments. We also comment on the 
case of a s-wave gap. 

Q , The nature of the carriers in the cuprates is an important question, which has been discussed and debated exten- 
^ ' sively. The availability of relevant experimental data over the years, and in particular angle-resolved photoemission 

(ARPES), should be helpful in the effort towards achieving a better understanding of this question. 
rS |- Here, we look at a characteristic one-particle property, relevant to ARPES. We obtain analytically the scattering 
^ rate in the superconducting state for a Fermi liquid with strong density of states peaks - van-Hove singularities (vHs) 
in 2-D - located close to the chemical potential fj,. Incidentally, this model does not apply to the underdoped cuprates, 
as there is no indication that a Fermi liquid description is appropriate. The qualitative nature of our results does 
not depend on the doping level, for as long as we remain within the realm of Fermi liquid theory. Precise numerical 
^ ' calculations can yield the quantitative dependence of the theory on the doping etc. A review of related work in the 
' frame of the so-called van-Hove scenario has been given in-'^ . The pinning of the vHs close to /i seems to be a plausible 
^ explanation for the common characteristics of a good many cuprates, whose van-Hove singularities are located between 
O 10-30 meV below the Fermi surface^, as ARPES experiments have shown. A review of some calculations yielding the 

1 pinning of the vHs close to /i appears in"^. The present work extends our results for the normal state, obtained in^'^, to 
the superconducting state. The analysis presented below can also be done in the frame of a weak-coupling BCS-type 

C\| . approach. We opt for the intermediate to strong-coupling Eliashberg approach, which is relevant for the cuprates. 
^ ' Qualitatively, i.e. as far as power laws etc. are concerned in terms of energy and temperature, the answers are the 
£^ , same for the two approaches. 

We consider a generic dispersion appropriate for the cuprates, of the type t,t',t" 
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£k — ~2t{cos kx + cos ky) — At' cos kx cos ky — 2t" (cos 2kx + cos 2ky) . (1) 



^ This dispersion generates strong vHs located at the points (±7r, 0), (0, ±7r). 

In the context of the Eliashberg formalism, the Green's function in the superconducting state is 

ienZ{k,en)+ek- fi + X{k,en) , , 

= 1^) ' 

Ch with e„ = (2n + 1)ttT and 

o . 

^' P{k,en)^[ienZ{k,en)]^ -[ek- fi + X{k,en)? -D\k,en) . (3) 



The renormalization Z{k, e„), the energy shift X(k, e„) and the (singlet) superconducting gap parameter D(k, e„) are 
given by the Eliashberg equations 

ie„[l - Z{k, e„)] = T 2^ ^ , (4) 



X{k,en)^Tj2 



V{k -p,en- en') [ep - M + X{p, en')] 

P(p,en') 



(5) 



and 
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V{q,LUn) generates the diagonal self-energy - c.f. below - and Vp{q,uj„) is the pairing potential. These potentials 
depend on the quasiparticle Green's function through screening. In general, e.g. for spin dependent interactions etc., 
these two potentials differ. 

The diagonal self-energy is given by 

S(fc, e„) = ie„[l - Z{k, e„)] + X{k, e„) = T V{k - p, e„ - e„0 G{p, e„.) . (7) 

Our derivation relies on this last equation, which is valid quite generally in the frame of a Baym-Kadanoff (BK) 
conserving approximation^, irrespectively of the specific Hamiltonian and approximation thereof^. 

The quasiparticle scattering rate is twice the imaginary part of T,{k,e) = Si(fc, e) + iT,2{k,e). It can easily be 
shown^ that ImT,{k, e) is given by the following formula for real energy at finite temperature : 



/oo 
— ImG^{k -q,e- w)ImV^{q,w) {coth(a;/2T) + tanh((e - co)/2T)} 
-oo 27r 



(8) 



We have 



ImG''{k,e) = B{k,e)Im{—^-^]+F{k,e)Re[—^-^} . (9) 

Here a = eZ{k, e),b=ek -H + X{k, e),D = D{k, e), B{k, e) = Re[ (^f^) and e) = Im[ (^^) . Note 

that a — b = G~^{k, e) in the normal state. 
We see that the factor 

p., f2{k,e)D,{k,e)-Mk,e)D2{k,e) 
^"''^ [h{k,e) + D,{k,e)f + [h{k,e) + D2{k,e)f' ^ 

with e) = Refik, e), /2(A;, e) = 7m/(fc, e), f{k, e) = eZ{k, e) + ek -n + X{k, e), and D{k, e) = Di{k, e)+iD2{k, e). 
This term is proportional to the gap function, and hence it gives a negligibly small contribution upon summation over 
the Brillouin zone for a d^2_y2 gap Di{k,e) ^ Dio{e)d{k), d{k) = coskx — cosky. So we can ommit the second term 
in the r.h.s. of eq. (9). This is not the case for a s-wave gap though, where this term should be kept. 
B{k, e) is given by 

o.. . _ /f (fe, e) + /|(fc, e) + Di(fc, e)/i(fc, e) + D^lfc, e)/2(fc, e) 
^ [h{k,e)+D^{k,e)]^ + [f2{k,e)+D2{k,eW ' ^ ' 

We can make the approximation 

^"^1 ^-j,}--TT5ie + fi-ek-^i{k,e)-Di{k,e)) . (12) 

( a — b — L> J 

This approximation was made in''''^. Setting ImG^{k,e) equal to a delta function is a reasonable approximation for 
this purpose, in view of the typical sharp spike feature of ImG^{k, e) shown in^ - c.f. fig. 1 therein. This is in the 
frame of numerical self-consistent FLEX calculations^. Further, numerically IniG^{k, e) is very small compared to the 
band energy for small couplings, and the difference of Im'E(k, e), as seen in our numerical calculation, for small and 
large coupling constants is mostly quantitative rather than qualitative. Moreover, the approximation of eq. (12) has 
also been used in the context of the electron-phonon problem in^ (c.f. the equation right before eq. (21.29) therein). 
Further, omitting terms of order [uP' 9^i?eS(A:, e)] and [w^ d^ReD{k,e)] we obtain 

i:§{k,e)c^~^B{k-q,e-Wk-g,e)IrnV'^{q,Wk-g,e){cot]i{wk-gj2T) + t&nh{{e - Wk-g,e)/2T)} , (13) 

Q 

where 

_ e-F/x - Cfc - Si(A:,e) - Di{k,e) 
l-dJl,{k,e)-dMk,e) ■ 
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We write 



oo ,,(2n+l).„N 2n+l 
n=0 ^ ' 



where Vq^\^) is the n— th derivative of ImV^{q,u> = 0) with respect to u>. Eq. (15) is true in the context of a BK 

approximation, and V{q, lu) may contain vertex corrections. E.g. for RPA and FLEX wc sec that the overall behavior 
of ImV closely follows the carrier susceptibility(-ies) Xo{q,'^), which we assume to be a regular function of w, as 

ImV{q,Lj) =^ImXoi{q,i^) \ei{q,uj)\'^ . (16) 

i 

ImXo is odd in w, while \e\^ is even, with e{q,w) the dielectric function. In the FLEX approximation, i =spin, charge, 
in RPA i =charge. In other BK approximations, e.g. c.f.^, eq. (15) is valid as well, with ImV{q,u) obeying an 
equation similar to cq. (16). 

The derivation given below for the linear in energy and temperature behavior of S2(fc,e;T) is generic. It relies 
essentially on a large coefficient for the linear in energy term of ImV - i.e. of ImXo etc. - and the presence of van-Hove 
singularities near the Fermi surface. The result holds regardless of the dimensionality of the system. However, it is 
important that a significant part of the spectral weight be included in the strong peaks of the density of states lying 
close to fi, a situation which is particularly favored in 2-D. 

For sufficiently small V^^"^(0), Vn > 1, we obtain 

l^V^J^iO) wk-g,e » 2.2^ (2n + l)\ ■ ^^^^ 

q q n=l ^ ^ 

This relation is valid for Wk-q,e < Cc, where the latter is the characteristic energy beyond which the infinite sum on 

the right becomes comparable to the Vq^^ term {cc can be estimated by equating the Ihs and rhs of eq.(17) with n = 1 
only). We constrain ourselves to the regime where eq. (17) is valid, as it is precisely this regime which gives the linear 
in max{T, e} dependence of S2 in the normal state**'^. 
Then we obtain 

I]2(/c,e) ~ -i^yg(i)(0)B(fc-g,e-u;fc_,,,) Wk-q,e {coth{wk-gj2T) + tanh((e - u;fe_g,,)/2T)} . (18) 
q 

The sum over q is dominated by the (extended) van Hove momenta qo in the vicinity of (±7r, 0) and (0, ±7r). 
Now notice that the following inequalities are satisfied: 

e + M - ek-q, > , (19) 
ReY.{k-qo,e) < . (20) 

The former is true for all momenta k in the Brillouin zone except (possibly) for /c ~ 0. The latter is seen to be true 
for e.g. < e < 4t from numerical FLEX calculations - c.f.^''^. Combining the two we get 

e + M - efc-9o - -ReS(A; - e) > . (21) 

Also, 

a,Si(fc,e) <0 , (22) 

and 

|9eSi(fc,e)| > |a,Z)i(fc,e)| , (23) 
for (5 < e < Co^'^, where (5 is a small negative energy and < Co = 0{t). Overall eqs. (19)-(23) yield 

Wk-qo,e > . (24) 
We obtain a linear in T scattering rate in the limit T > e, 'Wk-q,e 
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S^(fc,e)~-T^yW(0)B(fc-9,e-w;fc_,,,) . (25) 

The coefficient of T is maximized when k — Qo hes along the diagonals of the Brillouin zone, and minimized close to 
the vH points. This effect is of order (Di/fi)'^ though, so it does not contribute much for a small gap. 

It has been known long ago - sec c.g.^ - that the scattering rate becomes linear in T for T > 0)^/4, with ujb being 
the characteristic boson frequency mediating the carrier interaction. Our treatment shows that lob here is nothing 
else but the fermionic energy Wk-q^.e- The results here complement our earlier normal state results in'^'^. 

ARPES experiments^"'^^ have indicated that in the superconducting state, and, presumably, for e ^ 0, the scattering 
rate is linear in T for kp close to the nodal directions, but otherwise saturates to a value which increases as kp 
approaches the vH points. Our result is in agreement with ARPES for kp along the diagonals of the Brillouin zone. 
The reason for the off-diagonal discrepancy is not clear at present. Certainly more data, and on compounds othen 
than BSCCO, would be helpful in clarifying this issue. 

Wc now turn to the limit T < e. We note that for e < Wk-q^e and for w^-q^e < the contributions of tanh and coth 
annihilate each other in the low T limit. Then, for e > Wk-q„,e 

where Ck,e = 1+ < tanh((e — Wk-q„,€)/2T) > is a factor between 1 and 2. This equation yields a linear in e scattering 
rate, if dfT,i{k — Qo, e) varies slowly with e. Actually, assuming that Df{k, e) ^ ffik, e) we obtain 

E,(fc,e)^-^_^^^^^^,^^^|e + M-efe,-Si(fc,e) ffik',u.) + fi{k',w) / " ^'^^ 

where k' = k — qo, w = e — Wk-q„,e and No stands for the q-summation around a single vHs. We see that this 
expression is linear in e, and that the gap-dependent term is minimized for k — qo lying along the diagonals (with the 
assumption that fi{k',w) vary more slowly with k' than the gap). For a s-wave gap we obtain a dependence e — D 
instead, just from the term Wk-q^^e in eq. (26), and B{k,e) yields an additional D contribution. These two distinct 
dependencies of S2 on the gap arc characteristic of the respective gap symmetries. Estimates of the parameters in 
eq. (27) cannot be easily obtained analytically, due to the complicated nature of the Eliashberg formalism. Precise 
numerical estimates are possible, but they are beyond the scope of the present work. 

ARPES experiments in^'^'^^'^^'^'^ indicate that the scattering rate is a linear function of energy both in the super- 
conducting and the normal state, in agreement with our results. The normal state result was obtained in^''*. 

Summarizing, within the Fermi liquid Eliashberg formalism, we calculate the one-particle scattering rate in the 
superconducting state of overdoped cuprates. We compare our findings with existing ARPES experiments. For 
T < e, the e-dependence obtained is consistent with expts. For T > e, the T-dependence obtained agrees with expts. 
for kp along the Brillouin zone diagonals. 

The author acknowledges useful correspondence with M.R. Norman on ARPES data. 
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